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Thu  fiiui  leport  u  an  account  of  the  work  completed  at  the  McDonnell  Douglas  Research  Laboratonet  on  Office  of 
Naval  Resraich  Contract  No  NOOOI4-7t>-C-04‘*4  under  the  following  three  phases: 

Phase  I  Viscous  Rowfields  and  Airframe  Forces  Induced  by  Two-Dimensional  Lift  Jets  in  Ground  Effect 
( I  February  1 976  to  28  February  1 977) 

Phase  2  Compressible  Viscous  Flowfields  and  Aiilrame  Forces  Induced  by  TwivDimensional  Lift  Jets  in  Ground 
Effect  ( I  February  1977  to  28  February  1978) 

Phase  3  Viscous  Flowfields  Induced  by  Three-Dimensional  Lift  Jets  in  Ground  Effect  ( I  February  1978  to 
28  February  1979) 

This  work  was  done  in  the  Flight  Sciences  Department,  managed  by  Dr.  R.J.  Hakkinen.  The  principal  investigator 
for  the  first  pha*  was  Dr.  D  R.  Kotansky.  the  principal  investigator  for  the  second  and  third  phases  wu 
Dr  W  W  Bower  The  flowfield  numerical  work  was  done  by  I>r  Bower.  Dr.  R.K.  Agarwal.  and  Dr.  C.R.  Peters.  The 
computer  graphics  risutines  were  written  by  Mr  P.W.  Bieideman  and  Mr.  G.P.  Tomaschke. 

The  program  nxmitor  was  Dr  R.E.  Whitehead.  Office  of  Naval  Research.  Arlington.  VA. 
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I.  INTRODUCTION 

Wiih  the  cufient  mleiesi  m  jei  *nd  fan-pt)w«icd  venic»l-ukeolf-and-Unding(VTOL)  miliury  aircraft,  there  has 
been  an  increasing  demand  fi>i  impioved  pertoiinance-piediciion  methods.  This  demand  is  greatest  for  techniijues 
to  predict  propulsion-induced  aerodynamic  effects  in  tlie  hover  mode  of  VTOL  Higlit. 

Thu  (ask  is  not  easy,  hoivevct,  since  the  hover  mode  of  a  VTOL  aircraft  is  characterized  by  complex  flow 
phenomena  Ambient  air  is  entrained  into  llie  litl  jets  and  the  wall  jets,  leading  to  an  induced  downflow  of  an 
around  the  aircratt  and  a  resulting  suckdown  force.  In  addition,  the  inward  jet  flows  collide  and  create  a  stagnation 
region  from  which  a  hot  gas  fountain  emerges  and  impinges  on  tlie  lower  fuselage  surface.  The  fountain  is  a  source 
of  positive  induced  forces  which,  to  some  extent,  counteract  the  large  suckdown  forces  near  the  ground.  However, 
the  fountain  also  heats  the  airframe  surface  and  can  result  in  the  reingestion  of  hot  gas  into  the  inlet. 

The  VTOL  ground-effect  flow  is  characterized  by  thice-dimensionality,  high  turbulence  levels,  compressibility, 
stiimg  pressure  gradients,  and  regions  of  stagnation-point  and  sepatated  flow.  In  order  to  gain  a  fundamental  under¬ 
standing  of  lift-jel  induced  flows,  including  cntiaininent  and  fountain  characteristics,  the  McDonnell  Douglas 
Research  Laboratories  fMDRl  )  has  over  the  jiasl  three  years  conducted  a  theoretical  investigation  of  two-  and 
three-dimensional  impinging  jets,  under  contract  to  the  Oflice  of  Naval  Research  (ONR). 

In  the  first  two  years  of  the  contract,  attention  was  ftKUsed  on  a  single,  planat,  unvectoied  lift  jet  in  ground 
effect,  illustrated  Khematically  in  Ligure  I .  The  jet  exits  from  a  slot  of  widtli  D  in  a  contoured  upper-surface  a 
distance  H  above  the  ground  plane  The  region  of  interest  extends  a  distance  W  on  each  side  of  the  jet  centerline. 

The  jet  flowfield  can  be  characteti/ed  by  the  following  three  regions  a  free-jel  region  in  which  the  flow  is  essentially 
the  same  as  that  of  a  jet  issuing  into  an  unbounded  region  of  the  same  fluid,  the  impingement  region  in  which  the 
flow  changes  direction  with  a  large  pressure  gradient,  and  the  wall-jei  region  in  svhich  the  flow  traverses  the  surface 
with  zero  pressure  gradient.  The  fluid  surrounding  the  jet  is  entrained  at  its  boundaries  in  all  three  regions.  This 
entrainment  results  in  an  induced  secondary  P  w  which  can  produce  suckdown  forces  on  the  airframe  undersurface 
and  an  associated  lilt  loss 
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liolateii  re^KHis  of  th«  flowficld  illusiiaied  in  Figure  I  can  be  analyzed  with  simphned  lechniquea,  fur  example, 
Ruhel^  has  predrcled  (he  ground-plane  pressure  distributions  in  (he  unpingement  region  using  an  incompressible, 
inviacid.  rotational  model  fur  two-dimensional  flow.  Althixigh  good  agreement  has  been  obtained  with  measured 
ground-plane  pressures,  this  approach  cannot  provide  the  fuselage  undersurface  pressure  variation  (hat  is  necessary 
to  calculate  the  suckdown  force. 

A  mure  rigorous  analysis  of  the  planar  impinging  jet,  based  on  solution  of  the  time-aveiaged  Navier-Stokes 
equations  in  conjunction  with  a  turbulence  model,  has  been  adopted  by  MDRl.  in  its  contract  work  with  ONR.  In 
the  First  year  of  this  effiMt,  an  isolated,  planar,  irKompressible  lift  jet  in  ground  effect  was  analyzed,  and  in  the 
second  year  (his  analysis  was  extended  to  compressible  flow  In  both  endeavors,  the  lime-aveiaged  conservation 
equations,  combined  with  an  appropriate  turbulence  model,  were  solved  in  stream-funcliun/viMticity  form.  A 
conformal  mapping  procedure  was  used  to  (lanslorni  the  nonrec (angular  physical  plane  into  a  rectangular  cumpu- 
lational  plane,  and  a  second-order-accurale,  rmite-diflerence  scheme  that  is  stable  for  all  Reynolds  numbers  was 
applied  to  the  transformed  differential  equations.  Tlie  discretized  system  was  solved  iteratively  using  point  relaxation 
for  a  range  of  the  dimensionless  parameters  which  characterue  the  problem 

In  (he  third  year  of  the  contract  work  with  ONR.  MDRL  extended  its  Navier-Stokes  analysis  of  impinging  jets 
to  three-dimensional  geometries  The  conOguialion  of  primary  interest,  shown  schematically  in  Figure  2.  conusts 
of  two  interacting,  initially  axisymmetric  jets  with  fountain  formation.  The  jets  exit  from  two  axisymmetric 
nozzles,  each  with  diameter  D,  whose  centerlines  are  separated  by  a  distance  L  and  wliose  exit  planes  are  a  distance 
l{  above  the  ground  The  vodth  of  the  legion  ol  interest  is  the  dimension  W,  measured  from  the  jet  centerline.  As  in 
(he  case  of  the  planar  jet,  the  flowl'ield  can  he  characterized  by  free-jei,  impingemeni,  and  wall-jei  regions.  In 
addition,  a  fountain-upwash  region  is  formed  as  a  result  of  the  interaction  of  the  opposing  wall  jets  The  inclination 
of  (he  fountain  is  dependent  on  the  geometry  of  (he  configuration  and  the  ratio  of  (he  jet  strengths.  Ttie  fountain 
u  a  key  element  of  (he  problem  and  the  major  motivation  for  analyzing  this  geometry,  for  as  Wallers  and  Henderson^ 
conclude  in  an  auessmeni  of  current  meihods  to  predict  the  aerodynamics  of  V'STOI  aircraft,  "a  major  deficiency 
exists  both  in  available  data  base  and  prediction  methods  relative  to  fountain-induced  effects  " 
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An  inciwprfuiWe.  inviKiJ.  lotational  flow  model  that  is  applicable  to  three-dimensional,  jet  impinifeineni 
ptohlemi  has  been  described  by  RubeP  This  technique  has  been  used  to  compute  in  a  cartesian  coordinate  system 
an  axisy  minetrK'  jet  impinging  perpendicular  to  the  ground,  but  calculations  of  interacting  jets  with  fountain 
formation  have  not  been  reported. 

The  basis  of  the  MDRl.  appnach  adopted  in  the  ONR  contract  for  predicting  three-dimensional,  impinging  jct 
flows  IS  again  a  rmite-diffeience  solulum  of  the  time-aveiaged,  Naviei-Stokes  equations  with  a  complemeniary 
turbulence  model  The  conservation  equations  were  solved  in  terms  of  scalar  and  vectoi  potential  functions  and 
voriiciiy ,  Since  the  solution  domain  is  a  parallelepiped,  a  conformal  mapping  was  not  required  SolutuHis  of  the 
govctnmg  equations  were  obtained  with  an  extension  to  three  dimensions  of  the  second-ordei-accuiate,  linite- 
dilleience  algorithm  used  in  the  planar-jet  analysis.  An  iterative  solution  of  the  discieti/ed  system  of  equations  was 
obtained  with  line  lelaxatuMi  for  selected  dimensionless  parameters  which  define  the  problem. 

fills  lejvirt  describes  the  theoietical  approach  that  has  been  adopted  m  tlie  ONR  contract  work  to  compute  two 
and  three  dinseiisional  impinging  jets  The  flowfield  model,  numerical  solution  scheme,  and  representative  results 
which  illustrate  the  computed  entrainment  and  fountain-upwash  chaiactciistics  are  presented  The  diiection  ol 
tutuie  MDRl  work  in  the  calculation  ol  ground-effect  flowfields  is  diwussed. 
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2.  ANALYSIS  OF  TWO-DIMENSIONAL  IMNNCING  JETS 

In  thi$  section,  the  pioceduie  is  Jesciibed  foi  cumputinit  the  inconipreuible  ind  cumpiessible  turbulent  fluwfields 
for  the  unvectoted  planai  lilt  jet  in  ground  effect  shown  in  Figure  1  Contour  plots  of  the  calculated  flow  variables 
and  comparisons  with  experimental  data  are  presented 

2. 1  The  Governing  Equations 

The  iinse-aveiaged  continuit> .  momentum  (Naviet  Stokes),  and  thermal  energy  equations  for  steady  compressible 
tlow  are  given  below  in  tensor  notation  for  a  space-fixed  reference  througli  which  the  fluid  flows. 

Conservation  of  mass 


dtpu  ) 
- ^  -0 


(1) 


Conservation  of  momentum  m  the  lUl  direction 


dtpu  d.  ♦  pu  ”  u  "l 
_ 1  ) 

d*. 


(2) 


Conservation  of  thermal  energy 


dtpu  J  ♦  PU|"h”)  j 

dx!  ”  dT 


pu,  Uj  )  t  u,  (r,j--pu,  u,  ) 


dx. 


(3) 


In  these  equations. p  denotes  density,  u  the  veUxity  component  m  the  lUl direction,  p  the  static  pressure,  r,j  the 
shear  stress  in  the  ifh  direction  on  a  surface  normal  to  the  jf!}  direction,  g,  the  acceleration  due  to  gravity  m  the 
ith  direction.  J  the  total  specific  enthalpy .  h  the  statx  enthalpy .  and  the  heat  flux  in  the  jtl?  direction 

Fquations  ( I )  through  ( 3)  are  written  m  the  form  derived  by  Rubesin  and  Rose^  in  whxh  p,  p,  r,j,  and  q.  are 
represented  using  the  usual  Reynolds  decompositum  and  u,,  J,  and  h  are  represented  using  the  mass-averaged 
variables  introduced  by  Favre-'  Specifically,  in  the  Reyntdds  decompcnition,  a  general  flow  variable  $  is  expressed 
as  the  sum  of  a  time-averaged  component  ^  and  a  fluctuating  component 

(4) 

For  the  case  of  compreuible  flow,  this  approach  is  not  the  best  for  all  variables  because  the  resulting  equations 
contain  nurrKrous  explicit  mean-product  fluctuation  terms  arising  from  compressibility  and  do  not  have  the  same 
form  as  then  incompressible  counterparts  If,  however,  certain  variables  are  written  as  mass-averaged  variables,  the 
resulting  equations  have  virtually  a  term-by-term  correspondence  with  the  constant-property  case  In  the  Favre 
system, 

♦  (5) 

where  the  tilde.  **’,  denotes  a  mass  average  defrned  in  terms  of  the  conventional  tune  average  and 

#  -  #  -  .  (6) 

P 
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The  time-averaged  cimservalion  tqualuHis  ( 1 )  Ihrougli  (3)  have  the  same  form  as  the  instantaneous  cimservition 
equations  with  the  eaireptioii  of  the  Rey  nolds  stress  term  pU|  Uj  ,  the  Reynolds  heat  tlux  term  puj  h  ,  and  the 
mean-energy  dissipation  tetni  U|  (fij  l  /2pU|'u|  )  Moreover,  these  equations  essentially  have  the  desired  term- 
by-term  correspvHidence  with  the  incompressible  lime-averaged  conservation  equations  This  correspondence  can  be 
signiticant,  as  pointed  out  in  Reference  4,  if  the  basic  turbulence  mechanisms  in  the  compressible  flow  reiiuin 
sinular  to  those  in  the  constant-proper ty  flow  (the  so-called  Morkovin  hypothesis!  such  that  turbulence  modeling 
infointation  can  be  taken  from  incompressible  flow  data  and  applied  to  the  compressible  case 

In  order  to  obtain  a  closed  system  of  equations,  the  Reynolds  stress  term,  the  Reynolds  heal  flu*  term,  and  the 
iiKan-enetgy  dissipation  term  must  be  evaluated.  In  initial  work  on  the  planar  jet,  reported  in  References  6  and  7, 
\Nie-equalion  models  were  used  to  represent  the  turbulence  field  (the  Wollshlein  model*^  for  the  incompressible 
(el  and  Rubesm's  ealension  of  the  Glustiko  nuHlel**  for  the  compressible  (el)  Both  o(  these  approaches  suffer 
because  the  turbulence  length-scale  distiibuiKins  must  he  specified  For  this  reason,  more  recent  analyses  of  tfie 
incompressible  and  compressible  planar  (els  have  used  the  iwo-equalion  Jones-fjundet  turbulence  modellO  in 
which  the  length-scale  variations  are  implicit  in  the  analysis  and  need  not  be  specified  a  pnorx 

V^Tien  Fqualions  ( 1)  and  (3)  are  expanded  from  tensor  notation  into  cartesian  coordinates  (x,y)  with  the 
corresponding  velocity  components  (u,v),  the  result  is  the  so-called  primitive-variable  form  of  the  time-averaged 
conservation  equatums  for  a  compressible  planar  flow  In  the  present  work  the  equations  are  not  solved,  however, 
in  terms  of  the  prinutive  variables  (velocity  components  and  static  pressure).  The  reason  for  this,  as  indicated  by 
Ghia  et  al  in  Reference  1 1,  is  that  the  pressure  does  not  appear  as  a  dominant  variable  in  the  system  of  equations 
Since  only  the  first  derivatives  of  pressure  svith  respect  to  each  of  the  cwrdinates  ate  nested  in  the  system  and  not 
the  second  derivatives,  it  is  difficult  in  primitive-variable  solutions  of  the  Navier-Stokes  equations  to  maintain  the 
elliptic  character  of  the  pressure  field 

('onsequently,  in  the  present  approach,  the  governing  equations  are  solved  in  terms  of  a  dimensionless  compressible 
stream  function  v  defined  by 


—  ■  pu 

dy 


(7) 
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and  a  dimensionless  voiticity  u;  defined  bv 
()v  i)u 

ui  -  —  — 

if*  dy 

With  this  approach  the  stream  function  is  computed  Irom  an  equatiim  derived  b\  combining  l-quations  ( 7)  through 
(d(.  and  the  vorticily  is  computed  from  a  transport  equation  derived  by  combining  derivatives  of  the  momentum 
equations  m  a  manner  that  eliminates  the  pressure  from  the  equation  for  ui 

The  dimensionless  governing  equations  in  the  present  planar-flow  analysis  arc  given  below 


Poisson  equation  for  stream  function 
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in  suhsrqueiK  equaiums)  air  delinrd  in  (hr  Appendix  Die  h<.Hindai>  coiidilions  on  hquadons  ( 101  and  (ID  lollow 
lioni  known  vriocils  disdihulions  al  (he  houiidaiies  ot  (hr  soludim  domain 
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where  k  is  the  (urhulenl  kinelic  enei|0 
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where  k^n  is  (he  riles  (ive  (heimal  condiKlivi(x  and  s'p  is  (he  cons(an('piexsu(e  specilK'  hea(  Dir  houndarx 
condilions  on  (he  s(a(K'  en(halpx  are  drirrmined  (hnnith  ihe  lrmpeta(urr  and  hear  flux  disrrihudons  sprcilird 
a(  (he  Niundaries 

Transpsxr  eipialion  hn  lurhulenr  kine(K  rneiKx 
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l-qujiloni  ( 10)  thtoutih  ( 22)  contjin  lltr  MIowing  pataiiKlrtt  wtiicli  appeal  at  a  icsult  uf  the  noimali/alion 
Re  •  PiS’oD  jJi).  tl>e  Re>iu)IJt  nuniKei .  Ki  ■  Vjy  v^jto  ■*.  <1’^  Fi»uJe  number,  and  Pi  ■  Cp.  lnf,_o.  <be  Piandll 
numbet  Die  leleience  conditums,  denoled  b>  the  tubwiipl,  o.  aie  ta)ten  lo  be  ibe  cenleilinc  values  at  Ihe  )et 
exit  plane 

Fill  the  cate  of  iiKompieiiiblc  flow,  p  ■  I .  and  ibe  pieviout  equation!  become  ftreatly  vimplitied  both  in  form  and 
number  The  tianspoit  equation  for  iheinial  encif^' .  tbe  equations  of  state,  and  the  temperature  variations  of  the 
lianspoit  piopertiet  are  not  needed  In  addition,  the  Poisson  equations  foi  static  piestuie  need  not  be  solved 
umultaneouslv  with  the  Poisson  equation  foi  stieam  function  and  Ihe  voiticily  lianspoit  equaluMi 

2.2  The  Numerical  Solution  Scheme 

To  solve  the  llowfield  equations  lot  a  planai  je'  dischaigmi!  from  a  curved  upper  surface,  as  shown  in  Figure  I.  a 
conloima!  mapping  technique  is  used  to  map  the  iiiegulai  physical  plane  into  a  teclangulai  computational  plane 
The  governing  equations  are  then  rewiilten  in  leims  ol  Ihe  computational  plane  ccKiidinales  and  solved  using  finite- 
dillerence  procedures  .Since  there  is  geometric  lymmeiry  with  respect  to  the  jet  centerline,  isnly  half  of  Ihe  flow- 
field  need  be  treated 

T)ie  mapping  lechmque.  which  it  an  inverse  pnscedure.  was  devised  at  MI)Rl  by  G  H  llolfmaii  and  it  best 
explained  with  reference  lo  Figure  i  Inilially .  a  liriile -dillerence  compulalional  plane  with  ccKirdinales  ((.  f))  is 
specified  The  distance  fselween  nodes  in  the  (  direction  it  a  and  in  the  rj  direction  is  b.  where  a  and  b  are  not 
necessarily  equal  Stretching  functions  can  then  be  applied  in  each  ciKsidinale  direction 


Computational  ptarsc 
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Physical  plarsa 


Pifurt  3.  Iitvaraa  confoemal  mappint  ««lth  itratehlnt- 
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*»  ■  fj  (l». 


(23) 


ji)J 

A -I;  (7)  (24) 

Wii)i  i)icK  rcUtums,  i  nuppiii);  planr  in.  X)  is  Jcteiminrd  w(iic)i  permits  finer  resolutiun  of  the  fluwf'ield  in  reunions 
vstieir  t)ie  itradienis  of  t)i(  computed  vaniblrs  are  seveie  Finally ,  a  conformal  mapping!  given  by 

*i  ♦  iX  -  Fix  ♦  ly)  (25) 


IS  introduced  «(iich  specifies  the  physical  plane  (x,  y  1  Since  the  rrupping  (Fquation  (25)|  holds,  it  tollows  that 
both  X  and  y  obex  Laplace's  equations  with  ciHiidiiiate  stretching 
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The  boundary  conditions  imposed  on  Tquations  ( 2b)  and  (27)  follow  from  physical  constraints  when  they  are  known 
at  the  boundaries,  and  from  integration  ol  the  Cauchy-Riemann  lelations. 
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dr)  3y 

d*i 
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dX  3r) 

and 

dr) 

dt 

dX 

dii 

when  the  x  and  y  boundary  distributions  are  not  known 


(2S) 


(2«) 


The  unique  feature  of  the  present  mapping  scheme  is  that  instead  of  specifying  the  coordinate  distributions  in 
the  phvsical  plane  and  accepting  whatever  computational  plane  results,  a  nurrseiicallv  convenient  computational 
plane  is  specified,  and  the  corresptsnding  coordinate  distributions  within  the  specified  phy  sical  bsnindaries  are 
computed  For  this  reason,  the  mapping  schertK  is  an  inverse  priwedure 

In  the  present  analysis  the  mapping  equations  are  discreti/ed  using  conventional  centtal-difTerencc,  finite- 
dirtereiKe  approximations  to  the  derivatives  This  diKretiaation  is  not  the  case  for  the  flsiwfield  equations  The 
reason  fsH  this  can  be  explained  with  reference  to  the  following  general  form  of  the  flowfield  equations  which 
appear  in  the  analysis,  where  ^  is  an  arbitrary  flow  variable.  ^  p.  cu.h,  k.  or  rj 

a  -  ♦  >  —  *3  -  *  b  —  *  a  (30) 

xr 

The  coefficients  a.  7, 3.  *bd  ft  and  the  source  term  a  are  (unctions  of  the  flow  variables 
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riic  ccniial'Uiltcirncc,  riniic-Jiftcicncc  ilfioiithm  applied  lo  hquation  (.)0)  piovide^  dia|(iHial  doniinamx  and 
(hcfcbv  convcif^cnl  soluiions  loi  ih)I>  a  limiled  lange  in  (he  nugniiudei  of  the  cocfficienii  7  and  6  Thii  ranite  it 
(tenetallt  exceeded  in  (he  (raiiipoK  e<{ua(iont  in  which  7  and  6  contain  (he  Ke>noldt  number  as  a  multiplicative 
(actor  To  ensure  duftoiial  dominance  for  all  values  of  7  and  6,  Hoffman's  auiunented-central -difference  algorithm  I  ^ 
w'as  applied  to  Tquation  (30)  The  central  idea  ot  this  method  can  be  explained  by  ctmsidering.  lor  example,  (he 
derivative  i3o<dt  I'sing  the  point  of- the-compass  notation  sliowii  in  Tigure  4  lo  designate  the  vaiiable  locations,  this 
derivative  can  be  evaluated  at  the  central  point  P  in  teims  of  the  adjacent  east  (t)  and  west  (W|  points  with  (he 
following  truncated  Taylor-senes  tepiesenlalion  and  standard  ceniial-diflerence  approximation  to  (he  first  derivative 


a’  uV 

a  dp 

p’ 

df' 

P 

.1'  d{^ 

fl 


fi9urt  4  Fivt  poini  limt*  diffartnea  stanal. 


In  the  augrnentedcential  ditleience  scheme,  the  denvalive  <)(■'  is  retained  and  expressed  in  terms  of  lower-order 
derivatives  by  differentiating  the  governing  transport  eijuation  with  respect  to  {  Tlie  convective  derivative  rtd  rlrj  is 
represented  in  an  analogous  lasliion  with  the  algorithm 

The  finite-difference  eijuations  are  solved  using  point  relaxation  For  each  flow  variable  d.  a  local  residual  is 
defined  by  ■  'CNv  I  ^  where  M  is  the  iteration  counter  Convergence  lasaid  tube  achieved  when  r^  <  10'^ 
everv  where  in  the  flowfield 

2.3  The  Computed  Flowfieldx 

The  configuration  conudered  fint  in  (he  analysts  is  the  jdanai  impinging  jet  illustrated  in  Figure  I .  However,  in 
order  toaviad  the  additional  complexity  of  treating  (he  irregular  boundary  of  the  umiersurface  from  which  the  jet 
discharges,  the  top  surface  of  the  computational  repon,  as  shown  in  the  shaded  portion  of  Figure  5,  was  initially 
taken  to  be  a  distance  h  abiwr  the  ground  plane  The  right  surface  was  taken  to  be  a  distance  w  from  the  jet  center- 
line  SiiK'e  normal  impingement  is  considered.  geometrK  symmetry  about  the  jet  centerline  exists  so  that  only  half 
the  flowfield  need  he  solved. 

The  required  boundary  conditions  for  this  problem  follow  from  known  velocity  and  turbulence  profilet  at  the  top 
Nnindary  of  the  computational  repon.  from  the  noslip.  impermeable  wall  constraint  at  the  ground  plaiM.  from  sym¬ 
metry  at  the  jet  centerline,  and  from  the  assumption  of  no  gradients  in  (he  flow  properties  at  the  nght  boundary. 
CalculatHuis  were  pei formed  for  the  turbulent  planar  jet  with  Reynolds  numbers  of  1 .46  x  10^  and  1 .27  x  10^  baaed 
on  the  width  of  the  jet  J  at  the  top  surface  of  (he  computational  repon  and  the  jet  centerline  properties  at  the  ume 
station.  The  dimensions  h  ■  fi/J ■  I  and  w  ■  */3  *2.5  were  chosen  A  finite-difference  mesh  was  used  with  SI  nodes 
m  the  X  direction  and  41  in  the  y  direction. 
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Figur*  S.  Definition  of  the  computational  region  for  the  free-upper  boundary  planar  impinging  iat. 


Computed  How  vjtuhicj  for  Ihw  coni'igutjiion  with  Rc  »  1  27  x  10-^  arc  illuitrated  in  the  contoiii  plots  of 
Figure  fi.  riie  clurjcterniic  features  of  the  planar  impinging  je*  clearly  shown.  The  stream  function  distribution 
demonstrates  the  t  iii.ulaiion  /one  in  the  top  region  of  the  flowfield  and  jet  turning  near  the  ground,  and  the 
vortK'ity  distribution  demonstrates  the  convection  of  u)  towaid  the  right  boundary  with  the  devcIopnKiil  of  a  wall 
let  along  the  lower  surface.  The  entrainment  and  wall  jet  characlenslics  arc  further  illustrated  by  the  contour  of  the 
x-componeni  of  vehKily,  and  the  rale  of  decay  of  jet  centerline  veliKity  is  illustrated  by  the  contour  of  the 
y -component  of  veltwily .  Tlie  static  pressure  contour  shows  the  strong  pressure  gradients  expeneiKed  in  both 
coordinate  directions  when  proceeding  low-aid  the  stagnation  point,  a  feature  which  precludes  the  use  of  conven¬ 
tional  biHindary  layer  theory  in  modeling  the  viscous  effects  of  an  impinging  jet  flow.  The  turbulence  field  is  illus¬ 
trated  by  the  contour  plots  of  the  turbulent  kinetic  energy ,  the  turbulent  dissipation,  and  the  length  scale  The  latter 
IS  computed  from  the  relation  L  *  k*  -^./<j.  which  further  reflects  the  relative  turbulence  levels  in  the  flowfield. 

The  corresponding  contour  plots  for  the  planar  jet  with  Rc  »  1  46  x  1(7*  arc  not  shown  since  they  demonstrate 
the  same  qualitative  behavior  as  those  given  for  Re  ■  1.27  x  10^  in  Figure  6.  Figure  7  contains  comparisons  of  the 
measured  and  computed  ground-plane  pressure  distributions  using  the  data  of  Schauer  and  Eustis**  for 
Re  «  1.27  X  10-'  and  the  data  of  (>ardon  and  Akfiral'^  for  Re  ■  1.46  x  10*. 

In  caiculaiing  the  flowfields  for  the  planar  impinging  jel  with  a  confining  upper  surface,  the  computational  region 
shown  in  Figure  R  was  used.  The  solution  domain  has  a  height  H  and  a  width  W,  where  both  dimensions  are  normal 
i/ed  by  the  jet  width. 
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Fifure  •.  Property  contours  for  tfte  incompressrhfe  frM-upptr  boundary  planar  iat  in  grouitd  affact 
(h  •  1.  »  -  2.B.  Ra  -  1.27  ■  10*). 


D  Measured,  Schauer  and  Euttit^^, 
Re  ■  1  27  K  10^ 

o  Measured.  Garden  and  Akfirat^^, 
Re  -  1  46)(  10^ 

—  Computed,  Re  ■  t  .27  *  10^ 

—  Computed.  Re  •  1.46  x  10^ 


aenarreM 

Figure  7.  Comparison  of  experimental  and  theoretical  ground  plane  pressure  distributions  for  the 
normal  impingement  of  the  free  upper  boundary  planar  fat  mrilh  a  fully  deeeloped  profile. 


oereeMt  • 


Figure  S.  Definition  of  the  computational  region  in  the  physical  plane  for  a 
simulated  fuselage  undersurface  in  ground  effect. 


h'oi  the  case  ol  incoinptessible  How,  the  full»wtn|t  houndao  cunditiuns  ate  imposed  the  vurticity,  s'ream  func- 
tiiHi,  tuihulent  kinetic  eiierfty ,  and  tuihulent  dissipation  profiles  are  specified  within  the  entering  jet.  On  the  solid 
boundaries  conditions  result  from  the  constraint  of  a  no-slip  impermeable  wall.  On  tlic  jet  centerline  symmetry  con- 
ditums  are  imposed,  and  at  the  iiglit  boundary  the  flow  properties  ate  taken  to  be  uniform  with  respect  to  the 
tiansformed  x  oHudinate  Pressure  gradients  at  the  boundaries  are  specified  using  bquationsf  131  and  (141. 

Flowfield  calculations  have  been  made  for  a  lange  in  the  values  of  II  and  the  values  of  the  Reynolds  number  based 
on  jet  centerline  properties  at  the  exit  plane.  Kepreseniative  solutions  aie  sliown  in  this  section  for  the  upper  surface 
in  close  ground  effect  (II  •  2.  W  «  .1,68)  and  for  Re  «  1  .(X)  \  10*.  1 .00  x  1(H,  and  1 .30  x  10^  In  the  first  case  a 
linite  dilfetence  nKsii  was  used  with  41  nodes  in  both  the  x  and  y  diiections,  and  in  the  latter  two  cases  a  mesli  was 
used  with  81  m>des  in  the  x  direction  and  41  in  the  y  diiection. 

Contours  of  the  computed  fluid  propeiiies  aie  sliown  in  Figures  'i  and  10  lor  Re  »  1,00  x  10-  and  Re  »  1  .lOx  105, 
respectively  The  llowfields  fi>r  these  two  conditions  are  qualitatively  sinidai.  and  their  significani  characterislics  are 
the  following  the  stream  function  plots  illusliate  the  turning  of  the  let  and  the  recirculation  pattern,  including 
entiainment  into  the  primary  jet  Tlie  vorticity  distributions  show  the  convection  ol  u)  to  the  right  and  the  develop- 
iiKnl  of  the  boundary  layers  on  the  upper  wall  (where  sepaiation  iKCUrs  at  (he  point  of  uj  =  Ol  and  on  the  lower  wall. 
The  entiainment  and  recirculation  /ones  are  more  clearly  illustrated  in  the  plots  ol  (he  x-  and  y  velocity  components. 
The  static  pressure  contours  show  an  increase  in  p  going  toward  the  stagnation  jvnnt  where  the  value  is  set  to  anity 
Ot  ccHitse,  the  greatest  dilleience  in  the  liigli-  and  low- Rey  nolds  number  solutions  sliown  in  Figures  d  and  10  is 
retlected  in  (he  turbulence  tieldv  For  Re  •  I  30x  10-'  the  turbulent  kinetic  energy  and  length  scale  values  are  much 
larger  in  the  recirculation  /one  than  they  are  lor  Re  *  I  00  x  10- 

Tlie  contour  plots  for  the  planar  jet  with  Re  •  1  (K)  x  Ur*  are  not  sliown  since  they  have  the  same  qualitative  char¬ 
acteristics  as  those  for  Re  *  I  30  x  lo'  A  comparison  of  the  computed  and  measured  pressure  variations  tor  the 
ground  plane  and  simulated  fuselage  uiideisuilace  is  sliown  in  Figure  1 1  lot  the  Reynolds  number  range  with  H  *  2 
and  W  «  1  68  Tlie  data,  which  were  obtained  in  a  test  program  described  in  Ref  16.  were  measured  with 
Re*  I  .10  X  lo'.  while  the  computed  values  correspond  to  Re*  I  0(1  x  10-.  Ke  *  I  00  x  lO^.jnJ  Re  =  I  .10  x  lo'. 
The  influence  of  Res  nolds  number  on  the  solution  results  from  the  role  it  plays  in  deleiinining  the  turbulent  entrain¬ 
ment  into  the  primary  jet  A  comparivm  of  the  computed  and  measured  jet  lentetlme  luibiilent-kiiieticenergy  and 
veliKity  distributions  is  illustrated  in  Figure  1  2  for  the  same  geometry  and  values  of  Re  The  agreement  is  good 
except  in  the  proximity  ol  the  w.ill  where  the  Jones-1  aunder  turbulence  model  given  by  Fqualionst  Ihl  •(  ISJ  fails 
lo  provide  the  observed  variation  of  k 

For  the  case  of  compressible  flow,  the  following  additional  boundary  conditions  are  imposed  lo  determine  the 
transport  equalii'n  for  the  thermal  energy  within  the  entering  Kt  the  enthalpy  profile  is  specified,  the  upper  surface 
rs  assumed  to  be  isothermal,  and  the  lower  surface  is  assumed  lo  be  adiabatic,  and  on  the  jet  centerline  and  at  the 
light  biHindary.  the  enthalpy  gradients  with  resfiecl  lo  {  are  set  lo  zero 

Figure  1 3  shows  the  various  contour  plots  calculated  for  the  planar  compressible  jet  in  ground  effect  with  H  *  2, 

W  *  3  68.  Re  •  I  (K)  X  I0-.  Pi  *  0.<v8.  and  Fr  *  I  63  x  10-  Tlie  siieani-lunctioii.  voiticiiy,  and  velocity  vomjvoneni 
profiles  demonstrate  the  same  basic  form  for  the  compressible  case  as  for  the  mcompicssible  case  Tlie  static  pressure 
profile  IS  also  of  the  same  form,  but  p  at  the  stagnation  jvoint  is  not  set  lo  unity  since  its  value  is  fixed  by  the  equation 
of  stale.  The  static  enthalpy  ctmlour  shows  the  decay  in  temperature  with  distance  from  (he  jet  core,  and  the  density 
distribution  reflects  the  variations  in  both  the  static  enthalpy  and  pressure,  with  p  increasing  toward  the  upjver  bound¬ 
ary  where  'he  lemjveraiure  is  lowest.  The  k  and  t  profiles  (the  latter  is  not  show-nl  have  the  same  qualitative  character¬ 
istics  as  for  the  case  ot  incompressible  flow,  which  is  expected  from  the  Moikovin  hy  pothesis 

Higher- Reynolds-number  calculations  using  the  JoiKS-launder  turbulence  model  were  not  earned  out  for  the  com- 
pressble  planar  Kt  because  of  grid-u/e  limitations  imposed  hy  the  storage  capacity  of  the  CIX'  ('YBFR  175  used  in 
(he  computations.  These  calculations  must  he  perfiHmed  on  a  compufei  which  has  a  larger  storage  capacity  and  will 
thereby  permit  accurate  computatrons  at  higher  Reynolds  numbers. 
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Figure  11.  Comparison  of  eKperimental  arsd  theoretical  pressure  distributions  for  the  incompressible  planar 
let  in  ground  effect  (M  •  2.  W  ■  3.681 
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Figure  12.  Comparison  of  experimental  and  theoretical  centerline  property  distribuboiM  for  the  incomprewible 
planar  iat  in  ground  affect  IH  ■  2.  W  ~  3.681. 
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3.  ANALYSIS  OF  THREE-DIMENSIONAL  IMPINGING  JETS 


In  ihi^  »<ciu>n  tlnr  ptiKCiiuie  ib  ilcscribcii  loi  coniputin|(  the  inciimpicssihlc  (uihulcni  tlot^rielJ  resullin(!  t'lom 
two  intetactiii^.  initulK  jxis>ninietric  jets  with  lountain  toimaiion  as  sluiwn  in  l-igure  2  Contour  plots  or  the 
eak'ulateJ  flow  va'iables  are  presented 

3.1  The  Ciuveming  E(|ualiuns 

The  time  averagcd  continuits  and  nioinentum  equations  which  form  tlie  siartinji  point  in  the  three-dtiiieiisioiial , 
steads  ,  incompressible  How  analssis  are  i-quations  ( I )  and  (21,  which  can  be  simpliTied  with  the  condition  ol 
cimstant  densits  Turbulence  closure  is  obtained  bs  modeling  the  Resnolds  stress  term  using  the  (ilushko-Kubesiii 
one-equation  nuHlel  tor  incompressible  (low  two-equation  turbulence  model  iscurrenlls  nut  used  in  the  Ihree- 
diiiKiisional  analssis  because  ol  storage  limitations  imposed  bs  the  (A  BT  R  1 75  computer 

tquations  ( I )  and  (  2)  ate  expanded  Itom  tensor  notation  into  cartesian  C(Hitdinates(x.  y ,/ 1  with  the  corresponding 
vehicity  components  I  u.  v.  w ),  but  the  equatnnis  ate  not  solved  m  piimilive-vatiable  (otm  Rather,  a  poletilial- 
tunction  vorticity  approach  first  pioposed  bs  \iegbesola  and  Buries  ^  ^  is  adopted  In  this  technique,  winch  is 
analogous  to  the  stieam-tunction  voilicils  lepiesentation  of  the  two-dimensional  problem,  the  pressure  is  elimiiiated 
Itom  the  momentum  equations  m  the  loimulalion  ol  a  voiticils  liaiis|)oil  equation  The  pressure  field  is  then 
computed  liom  a  Poisson  equation  derived  Iroin  the  derivatives  ot  the  component  momentum  equations 

To  idenlicalls  satisis  the  conliiuiits  condition,  a  scalar  potential  ‘h  and  a  vector  potential  are  introduced  and 
the  velivits  field  is  given  as 


u,  *  —  • 

0X|  ■1'^  .-rx, 

wfieie  (||J^  IS  the  alternating  tensor  Hnasaki  and  Heliums  I  ^  have  shown  that  a  unique  vector  potential  \\.  within 
the  gradient  ot  an  aibitrais  haimonic  tuiKtion,  exists  such  that 


-\  ditlerential  equation  toi  the  vector  jsotential  is  derived  bs  taking  the  curl  ol  Tquation  (.52)  and  using  i  quatioii 
( -VH.  which  gives 


'  Wi  . 

.>x,.5x, 


where  the  vorticilv  is  detined  bs 


“■  ‘  .5x, 


A  differential  equation  for  the  scalar  potential  it  derived  bs  taking  the  divergence  of  Tquation  (.52). 
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Thr  jJvanijtic  ol'  iiiiiixlucint!  3  wjl;ir  potrniul  o  i)ir  convenience  ii  otTers  in  vetting  tlie  houniiar>  conJitiuns  on  the 
vriocilv  componenti  Since  hi^ujlion  (Aft)  is  nnlepciiiJeni  ol  KeynoIJs  iuinil>er,  it  need  be  solved  only  once  lor  a 
given  geometry  and  set  ol  boundary  conditions 

The  boundary  values  ol  the  vorlicity  are  evaluated  trom  I  quation  |35).  where  u,  is  calculated  liom  Equation  (32). 
The  bouiidaiy  conditions  on  A|  are  chosen  such  that  its  tangential  coniponents  aie  equal  to  zero,  the  conditions  on 
the  normal  comporieiKs  are  then  deterrmned  by  bquaitori  (33)  hot  a  Dal  boundary  the  conditions  become 

-  -O.A, -0.  (37) 

dll  ‘ 

where  n  denotes  the  normal  component  and  i  the  tangential  The  normal  velocity  components  aie  assumed  to  be 
known,  and  these  values  are  used  to  compute  the  normal  derivatives  ol  the  scalar  potential  on  the  boundary  , 

<>4> 

—  ■  nu  (38) 

dn  '  ' 


The  dimensusnless  governing  equations  in  the  present  iliiee-dimensional  Dow  analysis  ate  given  below 
Poisson  equation  lor  scalar  potential 

d'*  d‘<b  d‘<l> 

-  ♦  -  ♦  -  •  0 

^ 

d»*  dy*  d/* 

Poisson  equation  lor  x  component  ol  vector  potential 

d*A  d*A.  d-A. 

_ X  ^  ^ 

s  s  s  ' 

dx*  dy*  d/* 

Poisson  equation  lor  y  component  of  vector  potential 


(A*)) 


(40) 


d*Ay  d*A^  d‘A^ 


U. 


(41) 


dx*  dy*  dz* 

Poisson  equation  (or  z  component  of  vector  potential 


d*A^  d‘A,  d‘A^ 

■>  *  1  s 

dx*  dy*  dz* 


W. 


(42) 
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Traiupoft  «quati>.>ii  lot  turbuirni  luiictK'  cnfio 

/d*k  a-k  a*k\  I 
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Tho  ctinxijntx  <i.  !,>.  C.  xnd  X  hjvc  thf  vjliifv  at  0  2,  1  10.  a  ‘*3.  and  0  -4.  rcvjvi: livciv  and  inr  Rcvrnddi  luinibft  ot 
the  iiirbiikiivc  i  ix  dcl'infd  bv 


Tbf  Iftijilh  walr  dixtnbutiim  I  tnuxi  be  xpeiilied  thdiughoul  the  flow  Held 

rbe  Reynold}  number  ii  baxed  on  (el  icnierline  ptopeiliex  at  the  entrance  plane  lo  ibe  xoluinm  domain 
.V2  The  Niimrru'al  Srilulinn  Scheme 

In  the  xolulton  of  the  flow-field  equalionx  for  two  mteiactmg.  miiially  axixvmmeiiK  (clx  impinging  peipendiculaily 
on  the  ground,  a  conformal  mapping  ix  not  required  xinte  the  computational  domain  hax  no  iiiegulai  tvHindaiiex 
The  governing  equaliimx  tan  then  be  wntien  in  either  one  of  the  two  formx  given  below 

Poixson  equation 

V  V 

a*d  a*^  a*^ 

-  ♦ - ♦  -  •  a,  (501 

>  X  X  I 

ax*  ay*  a/* 


V 


Tiaiupott  equation 


d'0  3*0  3*0 

-  ♦  -  t  - 

■<  ^  ^ 

3a*  3y*  3/* 

wtieie  0  denotM  a  (teneial  floia  vanaMe  and  a.  d.  7.  *.  oi .  ay,  and  03  denote  gencttl  coelficient* 


\  do  3o 

)  ♦  Re  (}  —  ♦  Re  7  — 

I  d%  dy 


do 

Re  6  —  ■  Re  o,  ♦  o-, 
3/  *  ^ 


(51) 


The  Poiasiin  equations  are  diKietued  using  the  conventional  cential  diffetence  algorithm,  but  fot  the  reason 
discussed  in  the  section  on  the  planat  jet,  the  tianspoit  equations  are  discietiMd  using  the  augmented-cential- 
ditlerence  algonihm  due  to  Mottmanl-^ 

The  coupled  s>sieni  of  equatuHis  is  solved  iterativel)  using  Gauss  Seidel  line  relaxation,  where  optimum  relaxation 
factors  were  determined  on  the  basis  of  numencal  expeiiinentainm  As  in  the  planar  How  analysis,  the  convergence 
cnierion  lor  ea^h  variable  0  IS  based  on  a  residual  r^  ■  lOS+l  0nI<  10^’.  vshete  N  is  the  iteration  counter 


.V.)  The  t'ompuleJ  Row  fields 

Incompressible  thm fields  have  been  caUulaied  lor  the  three  dimensional  jet  configuration  of  Figure  2  using  the 
computational  domain  shown  in  I  igure  14  17ie  dimenu.uis  of  the  solution  domain  are  normalized  by  the  jet 
diameter  d  at  the  entrance  plane  to  the  computational  legion  Tlie  lop  surface  ol  the  latter  is  taken  to  he  at  a 
distance  M  h  below  the  ki  exit  plane  rather  than  at  the  exit  plane  in  order  to  decrease  the  size  o(  the  computational 
region  Three  ol  the  side  suilaces  are  jet  svmiuetrs  planes,  the  additional  svmineliy  plane  cut  through  the  fountain 
(z  ■  Ol  IS  not  made  so  that  jets  with  unequal  thrust  levels  can  be  treated  later  b>  the  same  computational  scheme 


Fifur*  14.  Oafimtson  of  tfs*  consputatsonal  rtfion  foe  two  intoeoclfn^  initiallv 
aaisymmatrsc  jttt  impenfinf  paepandscutorly  on  tfM  fround  iNith 
fountain  foemation. 


Brvundary  conditiisns  for  this  problem  are  determined  from  the  mvslip.  impermeahle-wrall  constraints,  symmetry 
conditiivns.  and  experimentally  determined  inflow  outflow  velocity  profiles  Turbulent-flow  cakulaticms  were  made 
with  I  «  .1.  h  »  I.  w  ■  I.  and  Re  ■  I  00  X  10*  In  setting  the  finite  difference  grid.  nodal  points  were  used  in  the  x 
direction,  f  in  the  y  direction,  and  25  m  the  z  direction  Use  of  the  TDC  7f»00  in  future  work  will  permit  a  greater 
number  of  mutes  and  iherehy  accurate  calculations  at  higher  Reynolds  numbers 


23 


Figure  I S  shows  rrprrseiilativr  Howfield  contours  in  various  planes  of  the  computational  region  sliown  in 
Figure  14  Plots  such  as  these  reflect  the  desired  symmelry  maintained  in  the  calculations  and  provide  the  basic 
data  requited  to  evaluate  entrainment  rates  and  fountain  pioperties.  Mowevei,  these  plana;  contour  plots,  even  in 
large  number  (which  ate  easily  generated  using  an  automated  plotting  priKedure  I,  do  not  provide  a  simple  visuali/a- 
Uon  of  the  overall  How  patterns 


The  physical  characterulics  of  the  llowfield  are  best  illustrated  through  thice-dimensional,  particle-pathhne  traces, 
as  shown  m  Figure  In  These  plots  are  constructed  in  the  following  manner  after  obtaining  a  convergent  solution 
of  the  Navier  Stokes  equations,  the  three  components  of  the  vcKx;ity  vector  are  stored  as  a  function  of  position  in 
the  computational  domain  A  point  is  then  selected  at  a  desired  pcMition  in  the  llowfield.  and  a  particle  is  moved 
along  the  local  velocity  vector  at  that  point  over  a  small  time  increment  When  the  particle  aiiives  at  a  new  position. 

It  IS  again  directed  along  the  local  veliK'ity  vector  This  priKedurc  is  repeated  stepwise  in  time  until  the  particle  exits 
the  computational  domain  Figure  lb  clearly  shows  the  particle  mosement  Irom  the  computational  domain  thtougli 
the  fountain  and  wall  jet  The  figure  demonstrates  the  capability  of  the  graphics  routine  to  rotate  the  flowl'ield.  which 
can  be  done  about  each  of  the  three  axes,  so  that  the  pathimes  are  visible  from  any  desired  angle 


ONR  has  funded  MORI  under  a  new  contract  \OtX)l4-7U<'.()h.t5  to  extend  the  flowtield  calculations  ol  ihiee- 
dimenuonal.  turbulent  impinging  jets  Additional  compulations  will  be  (seifoimed  for  equal-strength  jets  scilh 
normal  impingement  m  addition  to  compulations  for  ihe  following  cases  unequal-sliengih  jets,  normal  impinge¬ 
ment,  equal-strength  jets,  inclined  impingement,  unequal-strength  lets.  inclined  impingement,  and  an  axisy  mmeinc 
Kt  in  cmssflow.  The  accuracy  of  the  c«>mputed  llowlields  for  interacting  jets  with  noimal  impingement  and  fountain 
formation  will  he  determined  in  part  througfi  comparisons  with  data  which  will  be  obtained  in  an  experimental 
program  funded  by  ONR  and  conducted  by  the  Aenxly  namics  Department  of  McDonnell  Aiicralt  Company  under 
contract  NOOOI4  7q-C.0l.10 
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FlfHf*  IB.  ReprtMntathre  eontoun  for  the  iiwompretaible,  interacting  tats  in  ground  aHact 
(L  -  3,  h  ■  1.  *r  •  1.  Ra  -  1.00  n  10*). 


4.  CONCLUSIONS 


NklXinnell  IXHifclas  Research  Ijiboraiories  has  developed  an  approach  for  predicting!  two  and  three-dimensional 
turbulent  imptnfpnit  jet  flows  based  on  a  finite-difference  solution  of  the  fully  elliptic,  time-aveiafied  ciwiservation 
equations  of  fluid  mechanics  in  conjunction  with  a  turbulence  mixlel 

In  the  finl  year  of  a  contracted  research  program  with  ONR.  MORI,  developed  a  method  for  the  compulation 
ol  planar,  incompressible,  turbulent  lift  jets  in  the  proximity  of  airframe  and  ground  surfaces  The  equations  are 
cast  in  terms  of  stream  function  and  vorticity  and  are  solved  using  an  algorithm  which  provides  convergent  soluiimis 
for  a  wide  range  in  Reynolds  number  In  the  second  year  of  the  program,  MDRL  generalized  the  planar  jet  analysis 
to  include  the  effects  of  compressibility.  C'omparisims  between  theory  and  experiment  show  reasonable  agreement, 
although  additivinal  work  is  necessary  to  improve  the  turbulence  mtuleling  and  a  computer  with  storage  capacity 
greater  than  the  t'YBl  R  175  is  required  to  obtain  accurate  solutions  for  higlici  test  Reynolds  numbers  A  basic  con¬ 
clusion  drawn  from  the  planar-jet  studies  is  that  contounng  base  area  away  trom  lilt  nozzle  exits  and  perhaps  the  use 
of  nozzle  exit  extensions  will  significantly  reduce  aircratt  undersiirlace  suck-down  loads. 

In  the  third  year  of  the  ONR  contract  work,  MORI  developed  a  ptevedure  for  the  calculation  ot  three-dimensional, 
incompressible,  turbulent  lift  je  s  m  gii'und  effect,  with  emphasis  on  a  configuration  comprised  of  two  interacting 
jets  with  fountain  formation  ihe  equations  are  solved  using  a  polential-function  vorticity  formulation  and  an 
augnKnled,  finite-difference  algonthm  to  ensure  convergence  at  higti  Reynolds  numbers  Tlie  computed  llowfield, 
displayed  using  a  computet  graphics  routine  that  constructs  particle  pathlines.  demonstrates  the  expected  qualitative 
behavior  ll<>wevei,  additional  computations  and  comparisons  with  data  are  requited  to  quantitatively  evaluate  the 
ptivedure 

Further  wink  will  be  carried  <Hit  by  MDRl  under  ONR  contras  t  to  improve  the  accuracy .  efficiency,  and  general¬ 
ity  of  the  present  approach  tor  solution  of  the  Navier-Stokes  equatitHis  for  impinging  jets  As  Wallers  and  Henderson- 
stale  in  their  report  "V/STOl  Aerodynamics  Technology  Assessment"  "I  rom  a  heuristic  pinnl  of  view,  compuienzed 
analysis  and  prediction  techniques  for  V  STOl  aircraft  performanc-e.  particularly  in  the  hovenng  mode,  must  even¬ 
tually  lake  into  account  the  viscosity  and  turbulence  effects  inherent  in  the  Naviei-.Siokes  equations  "  WXile  Ihe 
NavKi-Siokes  analysis  developed  by  MDRl  cannot  provide  the  flowfield  associated  with  a  complete  VTOI  aircraft 
in  ground  elTecl,  it  can  provide  useful  information  on  the  jet  entrainment  and  fountain  Ilow  requiied  in  potential 
flow  paneling  methods  that  are  used  to  represent  complete  aircrafi  configurations  Data  of  this  type  cannot  always 
be  obtained  precisely  through  experiment  Moreover,  the  solution  scheme  described  in  this  paper  can  he  applied  to 
a  vanety  of  other  problerrw  of  interest  in  subsonic  aennly  namics.  such  as  flow  in  internal  passages. 
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APPENDIX: 

DERNITIONS  OF  FUNCTIONS  APPEARING  IN  THE 
GOVERNING  EQUATIONS  FOR  PLANAR  FLOW 

1  d\ii  I  dp  d~\li  dp  1  d^  I  dp  3*^  dp  3“|J/^ 

dx  \3x  dxSy  3y  3^2/  ^2  dy  \dy  3x3y  3x 

I  Id'rit  3-v;,  1  dp  ac-  ^  I  ap  d<;\  ^  z 

\  ax*  dxdy  I  ^3^2  3^2  p  dy  dy  p  d%  dxj  p  dxdy 


I  ap  dy 
p  dy  dx 


dkii  ^  t  dp  dii/  ^  d“C'\  -  ^^e(f  (/  Z  d^  dp  a*p\  dil/  dp  d*<i(  dp  3**1/ 

ax  p  dx  dy  ‘  dxdyj  -2  dx  \  y  p  dx  dy  dxdy/  dy  dy  dxdy  ~  ^  ^ 
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